Experimental evidence of paired hole states in model high- Tc compounds 
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The distribution of holes in Sri4_a;Caa;Cu2404i (SCCO) is revisited with semi-emperical reanalysis 
of the x-ray absorption (XAS) data and exact-diagonalized cluster calculations. A new interpretation 
of the XAS data leads to much larger ladder hole densities than previously suggested. These new 
hole densities lead to a simple interpretation of the hole crystal (HC) recently reported with 1/3 
and 1/5 wave vectors along the ladder. Our interpretation is consistent with paired holes in the 
rung of the ladders. Exact diagonalization results for a minimal model of the doped ladders suggest 
that the stabilization of spin structures consisting of 4 spins in a square plaquette as a result of 
resonance valence bond (RVB) physics suppresses the hole crystal with a 1/4 wave vector. 

PACS numbers: 



Two years after the discovery of high-Tc cuprates [l| , 
a new phase of Cu-0 based systems, Sri4_:rCaa:Cu2404i 
(SCCO), was found 0, Q. SCCO is a layered material 
with two different types of copper oxide sheets - one with 
Cu02 'chains' and one with CU2O3 'ladders' (see Fig. 
la). The sheets are separated by Sr/Ca atoms, and are 
stacked in an alternating fashion along the b crystallo- 
graphic direction. They are structurally incommensu- 
rate; the ratio of chain and ladder lattice parameters is 
Cc/cl«7/10. From charge neutrality, the formal valence 
of Cu is -1-2.25, resulting in 6 holes per unit cell if count- 
ing from a Cu^+ state. The substitution of Ca for Sr re- 
distributes the holes between the chains and the ladders 
allowing for studies as a function of hole density without 
the influence of strong scattering by charged impurities. 

A hole-doped two-leg spin ladder is the minimum 
needed to obtain superconductivity (SC) H, IM Qi al- 
though this competes with an insulating 'hole crystal' 
(HC) phase UHll- SC with Tc ~12 K has indeed been 
found by Uchara et. al Q for samples with x = 13.6 
under hydrostatic pressure > 3 GPa. Recently, the HC 
phase in the ladders was also discovered by resonant soft 
x-ay scattering (RSXS) [l3,|l3- The HC is observed to 
have only odd periodicity {\hc = 5cl for a; = and 
\hc ~ 3cl for X = 10, 11 and 12). The competition be- 
tween the HC and SC phases in the two-leg spin ladder is 
believed to be similar to that between ordered stripes and 
SC in doped two-dimensional Mott insulators [l^ mak- 
ing it an important model system to try to understand. 
In fact, a two-dimensional model of coupled two-leg spin 
ladders 13] was used to explain recent neutron scattering 
data of two-dimensional cuprates Lai5/8Bai/8Cu04 |14J . 

A serious problem exists however in understanding the 
periodicity of the HC in the ladders in terms of what 
was thought to be the correct doped hole density of the 



ladders. The often accepted doped hole density of 5l ~ 
1/14 (i.e ~ 1 hole for every 14 Cu's) in the ladder and 
5c ~ 5/10 (i.e. ~ 5 holes for every 10 Cu's) in the chain is 
a much too small density to arrive at a HC periodicity of 
bcL in ladder units. Until we understand this we really 
cannot reach any conclusions with regard to proposed 
models predictions regarding charge and spin ordering in 
these structures. 

The present paper deals with this problem and presents 
a new interpretation of the polarization dependent x-ray 
absorption data which leads to a much different charge 
density distribution between ladders and chains and with 
which the observed periodicities of hole crystals emerges 
quite naturally. We conclude that 5l and 5c are 2.8/14 
and 3.2/10 for x = 0, respectively, and the number of 
holes in the ladder increases almost linearly with x, reach- 
ing 5l = 4.4/14 and 5c = 1.6/10 for x = 11. The end 
result is a strong support of the rung based hole pairing 
predicted using simple t — J like models Q . We attribute 
the absence of the HC with A = 4cl, expected at a; = 4 
but not observed, to simple RVB-physics. 

The distribution of holes between the chains and lad- 
ders is essential in determining the answer for the prob- 
lem we mentioned above, and most physical properties. 
Like in high-Tc cuprates, the holes are expected to en- 
ter into O 2p orbitals and to form spin compensated lo- 
cal bound states with a Cu 3d hole il5| refered to as 
Zhang- Rice (ZR) singlets [13. Various experiments sug- 
gest that 5]^ « 1/14 and 5c « 5/10 at a; = 0, and provide 
evidence for a transfer of holes from the chains to the 
ladders upon Ca substitution. X-ray absorption mea- 
surements [131 find that 5h ranges from 0.8/14 for a; = 
to 1.1/14 for x = 12. Optical conductivity data [ig find 
a range from 1/14 to 2.8/14, while "^^Cu NMR studies 
^9^| find a range from 1/14 to 3.5/14. However, the opti- 
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FIG. 1: (a) Structure of chains and ladders. The orientation 
of the O 2p orbitals involved in the ZR singlets are indicated. 
The three different oxygen sites are identified for the ladder. 
The right-side panels show XAS for E||c and E||o and its the- 
oretical fits, for a; = (b,c), x = 4 (d,e), and a; = 11 (f,g). 
Squares are the experimental data. Black, red, green, cyan 
and blue lines, are theoretical curve for total fitting, respec- 
tively contributions of 0(ch), 0(1), 0(r) and UHB. 



cal conductivity and the NMR data analysis at cc = are 
based on the neutron diffraction observation ^O] of what 
was thought to be a superlattice reflection corresponding 
to a chain charge density wave (CDW) consistent with 
6c = 5/10, and thus Sl — 1/14. More recent neutron 
studies IJ, 12^ and analysis of the crystal structure by 
van Smaalen [23| have clearly shown that this peak is ex- 
pected in the basic crystal structure and is therefore not 
evidence for a CDW. 

The only direct measurement of the hole density dis- 
tribution comes from polarization dependent XAS. This 
is also subject to interpretation, and the model used pre- 
viously had unexplained discrepancies with regard to the 
polarization-dependence. In Niicker et. al's analysis |l7| 
of the XAS data it is concluded that the holes are mainly 
concentrated on the chains. Their interpretation assumes 
that there are only 2 distinct O Is pre-edge absorption 
energies, one (HI) corresponding to holes in the chains 
and the other (112) to holes in the ladders. The HI peak 
should be independent of the ac-plane polarization since 
the lobes of the O 2p orbitals involved in the chain ZR 
singlets are oriented at 45° to the a- and c-axes (see Fig. 
la). The H2 peak should be strongly polarization depen- 
dent since the rung O and the leg O have different hole 
amplitudes. However, for x — 0, XAS data shows that 
HI is as strongly polarized as H2. In Ref. 17] it is argued 
that this effect is small and therefore is neglected. 

Single crystals of SCCO were grown by travelling sol- 
vent floating zone techniques [2J| . The surfaces were pre- 
pared in the manner described in Ref. 10] . Polarization- 
dependent XAS measurements in the fluorescence detec- 
tion mode were carried out on the soft x-ray undulator 
line XIB at the National Synchrotron Light Source. The 
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FIG. 2: (a) Spectral weight (SW) (b) full width at half max- 
imum (FWHM) and (c) energy of 0(ch) (black), 0(l)c (red), 
0(1) a (green), 0(r) (blue), UHB for E||c (cyan) and UHB for 
E||a (magenta) as function of x. Blue open-circles are for 
oxygen ions in the ladders. 



energy resolution in the range of interest was about 200 
meV. The spectra were corrected for incident flux varia- 
tions and were normalized at about 70 eV above and 10 
eV below the edge where the absorption is atomic like 
and structureless. 

The a; = spectrum shown in Fig. 1(b) &(c) is identi- 
cal to that published in Ref. [13. Like them, we assign 
the lowest energy structures to the holes doped in O 2p 
orbitals and the higher energy structure at about 530 eV 
to transitions to the upper Hubbard band (UHB) or Cu 
3d orbitals. These are followed by a broader structure 
due to transitions to unoccupied bands hybridized with 
O 2p and 3p states. The UHB structure is only weakly 
polarization dependent, as expected given the symme- 
try of the empty (1^2 _y2 orbital. Since the point group 
symmetry for the ladder is not quite 04^ (the four O 
surrounding a Cu are not identical, see Fig. la) some 
polarization-dependence remains. In addition, there is 
a strongly polarization dependent feature at lower ener- 
gies, composed of at least two components. In Ref. 17], 
XAS data for La3SriiCu2404i.o2 (with only 3 holes per 
unit cell) shows only one nearly polarization independent 
structure at the lower energy. As concluded there, this 
strongly suggests that the holes involved are almost solely 
on the chains, where all O sites 0(ch) are identical, and 
with 2p orbitals oriented to 45° (see Fig. la). 

In ladders, things are more complicated. There are two 
types of O sites: the rung sites, 0(r), coordinated by 2 Cu 
ions, and the leg sites, 0(1), coordinated by 3 Cu ions (see 
Fig. 1(a)). The different coordination numbers result in 
different binding energies for the core Is and valence 2p 
orbitals. Higher values are expected for the orbitals of 
0(1), while those of 0(r) should be close to 0(ch) which 
is also coordinated by 2 Cu ions. Moreover, each ladder 
ZR singlet involves one 0(r)a, two 0(l)c from the leg, and 
one 0(l)a from the leg of a neighboring ladder. The sub- 
scripts a and c refer to the polarization needed to observe 
O Is to 2p transitions. For a-polarization, transitions are 
possible for 0(r)a and 0(l)a at different energies, while 
for c-polarization, transitions are only possible from 2 
identical 0(l)c, with energy close to that of 0(l)a. 

We performed a simultaneous least-square fit to all the 
measured spectra. Outputs of the fits, namely the spec- 
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tral weight (SW), full width at half maximum (FWHM) 
and various energies are shown in Fig. 2(a)-(c). The 
number of holes Sl — (3 ~ 5Sc)/7 are determined from 
the spectral weights of the various absorption lines: 

" 7{SWoir)^ + 5W^o(o„ + SWon)^ + SWoich)) ^ ' 

where SWo(ch) is the total SW for both polarizations. 

We start with x = 0. XAS data and theoretical fits 
are shown in Figs. 1(b), 1(c). For E||c, the doped hole 
region has at least two structures. In our interpretation 
this is due to the energy difference between 0(ch) and 
0(l)c. For E||a, the contribution of 0(l)c should be re- 
placed with that of 0(l)a and 0(r)£,, thus shifting more 
weight to the lower energy. This is indeed consistent with 
the data. In our fitting results, shown in Fig. 2, the en- 
ergy of 0(r)a is about 0.2 eV higher than that of 0(ch), 
while the energy of 0(1) c and 0(1) a are roughly equal 
and about 0.5 eV higher than that of 0(ch). The SW of 
the UHB is almost polarization independent. This shows 
that the ladder holes are distributed nearly isotropically 
amongst the 4 O involved in the ZR singlet, even though 
the symmetry is not the full D4/J. (If large deviations 
were found, the ZR picture would not be valid for the 
ladder holes). The FWHM of 0(ch), which is about 0.62 
eV, is about 5% larger than the ones in the ladders. From 
Eq. (1), we find Sl = 2.8/14 and 5c = 3.2/10. 

We continue the analysis for x — A. In Figs. 1(d), 1(e) 
we show the fitting for E||c and E||a data, respectively. 
Again, the SW in the UHB is almost polarization in- 
dependent (Fig. 2(a)) while the energies of the various 
O sites are close to the x = values (Fig. 2(c)), con- 
sistent with our understanding of the polarization de- 
pendent XAS. From Eq. (1) we find Sl = 3.4/14 and 
Sc = 2.6/10. 

Figs. 1(f) and 1(g) show our XAS data for x = 11. 
(The maximum HC intensity occurs at x — 11 where 
the wave vector is closest to 3cl [13 )• ^^ R,ef- 01 it 
is claimed that here, the UHB is strongly polarization 
dependent while the O doped hole pre-edge region is al- 
most polarization independent. Our results show an op- 
posite behavior, like for samples with smaller x. In fact, 
the shape and the intensity of the hole-doped and UHB 
peaks of Ref. [13 are similar to ours, and can be made to 
coincide by rescaling. We conclude that the difference is 
due to the method used in the normalization of the data. 

The fit of the polarization dependent XAS for x — 11 is 
more difficult than for cc = or 4. The problem is not sta- 
tistical noise, but rather the structural peaks themselves. 
We see in Fig. 1(f) & (g) that there is no clear evidence 
for multiple peaks in the doped hole regime because the 
energy resolution is too poor to resolve the peaks. We 
therefore use the a; = results as input for the fitting. 
The fits are shown in Figs. 1(f), 1(g). As before, we find 
that the SW of UHB is almost polarization independent 
and the energies of the various O sites remain close to the 
X — values (see Fig 2), validating our interpretation. 
We find Sl = 4.4/14 and Sc = 1.6/10. It is important to 
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FIG. 3: (a-i)-(a-iii) Scenarios for the ladder hole distribution 
in a HC with A = 4cl, as predicted by DMRG. Filled (open) 
circles represent electrons (holes) . Arrows indicate the spin of 
unpaired electrons, (b) The estimated number of holes in the 
ladders (filled circles) and the chain (open circles). Dashed 
lines show the number of holes for scenario (a-i), for Xhc/cl = 
3, 4, and 5. (c-i) Disordered and (c-ii) ordered hole arrange- 
ment corresponding to \hc = ^cl- The disordered state has 
equal numbers of 2-rung and 4-rung spin plaquettes. (d) GS 
energy per rung, in units of Jy , as a function of the number 
of rungs in the spin plaquette. Black, red, green, blue, cyan, 
and magenta-filled circles are results for J±/J\\ ~ 1.2, 1.0, 0.8, 
0.6, 0.4 and 0.2, respectively, (e) Edis — Sord as a function of 
Jx/J||, for J|| — 130 meV (see text for details), (f) easiNR) 
vs. Nr for J±/J\\ = 0.4 and Jri„g /J\\ = -0.1, 0, and 0.1 
(black, red and green-filled circles, respectively). 



mention that these results depend strongly on the energy 
of 0(ch). For example, varying it by 0.05 eV changes Sl 
by about 0.5/14. This is due to an instability in the fit- 
ting process caused by the close proximity of peaks in 
the pre-edge region (the energy of 0(ch) falls close to the 
leading edge of the hole-doped peak). This is why the 
X — spectrum provides an important reference. 

We now analyze some possible scenarios of the hole dis- 
tribution in the ladder, for these new ul values. We also 
consider the connection to the wavelengths Xhc = 3cl 
{x = 11) and 5cl {x = 0) of the recently discovered 
HC l^lnj. First, DMRG calculations for a single lad- 
der [3 found that holes prefer to pair along the rungs, 
resulting in a charge density wave shown pictorially in 
Fig. 3 (a-i). For periodicity Xhc = Ncl, the doped hole 
density in ladder Sl should be 1/N, i.e. Sl — 0.2 if iV = 5, 
Sl = 0.25 if iV = 4 and Sl = 0.333 if TV = 3. These vah 
ues are very close to our new XAS estimates, as shown 
in Fig. 3(b). Surprisingly, however, the Xhc = 4cl 
HC, which is predicted by DMRG [ 1 11| to be a stable 
phase, is not observed in RSXS ill|. While the study 
of this discrepancy is on going, a possible explanation 
is proposed below. Another hole distribution consistent 
with these Xhc values, not involving pairing, and also 
discussed in Ref. 0, is single rung bond-centered holes 
which is shown in Fig. 3(a-ii). Since Sl for this case is 
halved, this model does not match our XAS results. In 



a third scenario proposed in Ref. [3|, the holes could be 
site centered, alternating between the two legs as shown 
in Fig. 3(a-iii). While this matches the Sl values, it re- 
quires an odd number of undoped rungs in the HC unit 
cell and is therefore inconsistent with the observed A^ = 3 
and N — 5 HC. We conclude that our results support the 
scenario of holes paired along the rungs. 

Assuming rung-paired holes, a minimal model of the 
doped ladder is an antiferromagnetic Heisen berg Hamil- 
tonian plus a cyclic four spin exchange term \2n . l27i | : 

Nr-1 Nr 

^ ~ "^11 / ^ ^a,n ■ ^a,n+l + J± / ^ ^l,n ' S2,n + ^ring 
"=i n=l 



Here, Nr ^ 
unit cell, a 



A^— 1 is the number of undoped rungs per HC 
= 1,2 indexes spins on the two legs, and J± 



and J|| are exchange couplings along the rung and leg, 
respectively. We assume no coupling between spins on 
opposite sides of rungs occupied by paired holes. Tiring = 

■^ring 2_/n=l [ 4 +^1," ■Sl,n-|-1 + S2,n •S2,n+1 +Sl,n •S2,n+1 + 
S2,n • Si_„+i + Si_„ • S2,n + Si_„+i • 82,™-!-! + 4{(Si_„ • 
S2,n)(Sl,,i+l-S2,„+l) + (Si^„-Si_„+l)(S2,n'S2,n+l) — (Sl,„■ 
S2,n+l)(S2,n ■ Si^„+i)}] is a four-spin cyclic exchange. 

We use exact diagonalization to find the ground state 
for various Nr values. The ground-state (GS) energy per 
rung, ecsiNn) = Egs{Nb.)/Nr, is shown in Fig. 3 (d) 
for various ratios of J±/Jn and Jring=0. The value of 
J±/ J\\ is not known accurately, but is believed to be be- 
tween 0.5 and 1.13 J^, 29, 30, 31, 32, 33]. An even/odd 
oscillation is observed for small Nr and J±/ J\\ < 1, fa- 
voring Nr = 2 and 4 {Xhc = 3, 5cl). The origin of this 
oscillation is simple. The limit J|| ^ J± corresponds to 
two AFM chains weakly coupled along the rungs. For 
even Nr, spins on each leg pair in a RVB-like state, and 
Egs is low. For odd Nr, each leg has an unpaired spin, 
significantly increasing Egs- In the limit J± ^ Jy, the 
GS consists of spin-singlets along the rungs and the par- 
ity of Nr is irrelevant. At large Nr, cgs converges to the 



bulk value. This even/odd oscillation provides a possible 
explanation for the absence of a HC with Nr — 3 {Xhc — 
Acl)- This HC costs an energy Cord = 6eG5(3) per two 
unit cells (see Fig. 3(c-ii)). A disordered phase, at the 
same doping, has equal numbers of Nr — 2 and Nr^ — 4 
plaquettes and an energy Cdis — 2eGS'(2) -I- 4eGS'(4) for 
the same length (see Fig. 3(c-i)). If edis < &ordi the HC 
phase is unstable. For Jy — 130 meV, we plot Cdis — ^ord 
in Fig. 3(e), showing that the disordered phase is ener- 
getically favorable, especially for lower values of J_l/J||. 

We have also studied the effect of Wring on Eqs ■ Such 
terms appear in 4*'' order perturbation expansions in the 
strong coupling limit of the Hubbard model [2g and are 
known to play an important role for Wigner crystals and 
''He solid. Typical results for eGs(^fi) are shown in 
Fig. 3(f), for Jj_/J|| =0.56 ^j and Jring/^|| = -0.1, 
0, and 0.1. Since the sign of the ring exchange and su- 
perexchange should be the same [2g, it follows that a 
large Jring suppresses the even-odd effect. For example, 
for J±/J\\ =0.56, edis — Sord increases from -20.8meV if 
Jring = 0, to -12.5 meV if Jring/J|| =0.1. We conclude 
that for reasonable values of J±, J\\ and Jring this sim- 
ple model offers a possible explanation for the absence 
of the Xhc ~ 4ci HC. An accurate determination of the 
exchange couplings is needed before the issue can be set- 
tled. 

In conclusion, we propose a new interpretation of po- 
larization dependent XAS for SCCO. Based on our anal- 
ysis combining the XAS and RSXS data, we find strong 
support for a pairing of holes in the rungs of the ladders. 
We also give a possible explanation for the absence of the 
HC with 1/4 periodicity in terms of RVB physics. 
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